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ABSTRACT 

Recently, Planck measured a value of the cosmic microwave background (CMB) optical depth 
due to electron scattering of t = 0.066 ± 0.016. Here we show that this low value leaves 
essentially no room for an early partial reionisation of the intergalactic medium (IGM) by 
high-redshift Population III (Pop III) stars, expected to have formed in low-mass minihaloes. 
We perform semi-analytic calculations of reionisation which include the contribution from 
Pop II stars in atomic cooling haloes, calibrated with high-redshift galaxy observations, and 
Pop III stars in minihaloes with feedback due to Lyman-Werner (LW) radiation and metal 
enrichment. We hnd that without LW feedback or prompt metal enrichment (and assuming a 
minihalo escape fraction of 0.5) the Pop III star formation efficiency cannot exceed ^ a few x 
10“'*, without violating the constraints set by Planck data. This excludes massive Pop III 
star formation in typical 10 ®Mq minihaloes. Including LW feedback and metal enrichment 
alleviates this tension, allowing large Pop III stars to form early on before they are quenched 
by feedback. We hnd that the total density of Pop III stars formed across cosmic time is 
< xo'*”® Mq Mpc“® and does not depend strongly on the feedback prescription adopted. 
Additionally, we perform a simple estimate of the possible impact on reionisation of X-rays 
produced by accretion onto black hole remnants of Pop III stars. We hnd that unless the 
accretion duty cycle is very low (< 0.01), this could lead to an optical depth inconsistent with 
Planck. 

Key words: stars:Population III-galaxies:high-redshift-cosmology:theory 


1 INTRODUCTION 


Understanding cosmic reionisation is currentl y one of the most 
exciti ng frontiers in astrophysical cosmology dLoeb & Furlanett9 
l2013h . Constraints from Lya forest measurements indicate that the 
IGM was completely ionised around z ~ 6 — 7 jpan et al.ll20Q6l: 
iMortlock et alj|2oTlI : ISchroeder et alj|2013l : iMcGreer et al. l2015h . 


however the exact ionisation history at higher redshifts remains un¬ 
clear. The optical depth due to electron scattering of the CMB, r, 
provides an important constraint, however because it is only one 
number (an integral constraint on the ionisation evolution), it is 
degenerate with different ionisation histories. Future data, includ¬ 
ing radio ob servations of 21cm emission from neutral hyd rogen 
in the IGM jFurlanetto et alj|2006l : IPritchard & LoeblbOl^ . will 
give a more detailed picture of the reionisation process. Comparing 
this data with theory will provide an important test of the standard 
model of cosmology and yield information about ionising sources 
which may be too faint to observe directly. 

Reionisation is thought to be primarily driven by UV pho¬ 
tons from stars. The first stars in the Universe are expected to 
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form fr om metal-free gas in 10® Mq dark matter “mini¬ 
haloes” (tHajm anet alj 19961 : iTegmark et al.ll9^ : lAbel et alj2002l : 


ISromm et a iJ l2002l) . These so-called Population III (Pop III) 
stars are predicted to be more efficient at pro ducing ionising ra¬ 
diation than metal-en r iched (Pop I/II) stars (iTumlinson & Shulll 
l200d : ISchaerej I2002L 1200 3h . This ionising efficiency would be 
further enhanced if Pop III stars form with a top-heavy initial 
mass function (IMF). Cosmological hydrodynamical simulations 
have been utilised to study the formation of Pop I II stars (e.g. 
Stacyet^ ]|2010l:IClarketalJl201ll: iGreif et alj|2012l: iHirano et all 


20141 : [Brommll2013l : Greili 2015h . but the IMF remains highly un¬ 


certain. 


Formation of Pop III stars in minihaloes requires efficient 
molecular cooling. As stars form, a background of Lyman- 
Wemer (LW) radiation builds up over cosmic time. Eventually 
this LW radiation is strong enough to photo-dissociate molecu- 
lar hydrogen, inhibi ti ng additional star formation in minihaloes 
jHaimane^^ 11 19971: iMachacek et alj 1200 ll: IWise & Abell l2007l: 


(Tjhea & NormaJ 20081 : Wolcott-Green et alj 1201 ll : Visbal et^hl 

2014h . At this point, only haloes with virial temperatures Tvir > 
10* K (hereafter “atomic cooling haloes”) can form stars through 
atomic hydrogen cooling. Pop III stars could form in these haloes 
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2 E. Visbal et al. 


if they contain pristine gas, but once they become metal-enriched 
from previous generations of stars, Pop I/II star formation occurs. 

Recently Planck has reported an improved measurement of 
the CMB electro n scattering optical depth, r = 0.066 ± 0.016 
( IXd^ et alj|2015t). This is low er than previous measurements from 
WMAP jKomatsu et al.ll20llh and leaves less room for an early 
partial reionisation. In this paper, we examine the constraints this 
new measu rement puts on the production of Pop III stars i n mini¬ 
haloes (see iHaiman & Brvan|[2006l : lOreif & BrommlllOO^ for re¬ 
lated studies based on the WMAP cosmological parameters). We 
utilise a semi-analytic model that includes the contribution from 
both atomic cooling haloes, calibrated with observations of a « 6 
galaxies, and minihaloes with a self-consistent treatment of LW 
feedback. Using a simple treatment, we also include the possible 
effects of metal enrichment in minihaloes from Pop III supernovae, 
which causes a transition to Pop II star formation. 

We compute the ionisation history and corresponding CMB 
optical depth for a variety of model parameterisations with and 
without LW feedback and metal enrichment and find that without 
LW feedback or metal enrichment, massive Pop III stars cannot 
form efficiently in minihaloes without violating the Planck con¬ 
straints. When LW feedback and metal enrichment are included, 
massive Pop III stars could form efficiently early on, but they 
are suppressed at lower redshifts, reducing the optical depth suffi¬ 
ciently to be consistent with Planck. We also find that, irrespective 
of the feedback prescription used, the total density of Pop III stars 
formed over all cosmic time cannot exceed « Mq Mpc“® 

without violating the Planck optical depth constraints. 

Note that iRobertson et ^ j2015l) recently performed an anal¬ 
ysis of reionisation and the Planck optical depth. Their study is 
empirical, focusing on the ionisation history implied by the ob¬ 
served UV luminosity function (LF) (with a modest extrapolation 
to fainter galaxies), while we address the imp lications for lower - 
mass minihaloes and Pop III stars. Similarly, iMitra et ^ l l2015ll 
use a semi-analytic model to show that reionisation from Pop II star 
formation alone (i.e. no Pop III) is consistent with the Planck opti¬ 
cal depth measurement and high-redshift quasar absorption spectra. 
They do not attempt the put upper limits on the amount of Pop III 
star formation permitted, which is the primary goal of this work. 


In addition to the UV photons from stars. X-rays from 
blac k hole accretion could potentially contribute to reionisa- 


tion 

Venkatesan et alj200ll:lMadau et al.ll2004l;lRicotti & Ostriked 

|2004 

; Volonteri & Gnedinll2009b. We perform a simple calculation 


to estimate how much gas the black hole remnants of massive Pop 
III stars could accrete with a high radiative efficiency without pro¬ 
ducing a r inconsistent with the Planck measurement. We find that 
for our fiducial model with LW feedback and moderate minihalo 
star formation efficiency, black hole remnants of Pop III stars can¬ 
not accrete at the Eddington limit with a duty cycle higher than 
~ 0.01. This suggests that either massive Pop III stars are un¬ 
common or there is s ome feedback mechanism which preven ts sus¬ 
tained accretion (e.e. lAIvarez et alj2009l : lTanaka et al.l20I2h . 

This paper is structured as follows. In § 2, we describe our 
reionisation model including the self-consistent prescription for 
LW feedback, the various model parameters and their chosen fidu¬ 
cial values, and our simple treatment of metal enrichment from 
Pop III supernovae. We present the results of this model in § 3. In 
§ 4, we perform a simple calculation to estimate how the growth 
of black hole remnants of massive Pop III stars would impact 
reionisation. Finally, we discuss our results and conclusions in 
§ 5. Throughout we assume a ACDM cosmology consistent with 


the latest constraints from Planck jAde et ai]l2014l) : Ua = 0.68, 
Um = 0.32, Ub = 0.049, h — 0.67, <js = 0.83, and ris = 0.96. 


2 REIONISATION MODEL 

Here we outline our semi-analytic model of reionisation. In § 2.1 
and § 2.2, we explain how we compute the ionisation history and 
describe our self-consistent treatment of LW feedback. In § 2.3 we 
discuss the physical parameters of the model and their fiducial val¬ 
ues. Finally, in § 2.4 we introduce our simple treatment of metal 
enrichment due to Pop III supemovae. 


2.1 Ionised filling factor 

We model the global reionisation process by c onsidering star for¬ 
mation in dark matter haloes, closely followine IHaiman & Holden 
l l2003b . An ionising efficiency and associated ionised volume are 
assigned to each dark matter halo and the total halo abundance 
is computed analytical ly with the Sheth-Tormen mass function 
dSheth & Tormenlll999l) . We assign different ionising efficiencies 
to minihaloes which we assume host Pop III stars and atomic coo¬ 
ing haloes hosting Pop II stars. We also assume that in regions of 
the IGM that have already been ionised, the increased leans mass of 
the photo-heated gas prevents star formation below a characteristic 
halo mass. Mi. The other important mass scales are the minimum 
minihalo mass and the atomic cooling mass, which are denoted by 
Mm and Ma (we discuss the fiducial values of these masses be¬ 
low). It follows that the total ionised filling factor as a function of 
redshift, Q{z), is given by 


Q{z) = 


p\,{z) / dz 

J OO 


C^-fcoll,i / /\ , 


(1 - Q(. )) X ) 


V{z',z), 


( 1 ) 


where Pb{z) is the mean cosmic baryon density. The ionising effi¬ 
ciency in minihaloes (i.e. number of ionising photons escaping into 
the IGM per baryon incorporated into a dark matter halo) is given 
by Cm = /*, m/esc,mt7ion,m, where the star formation efficiency, 
/*,m, IS the fraction of baryons in mimhaloes that form stars, ^esc,m 
is the ionising photon escape fraction, and ?7ion,m is the number 
of ionising photons produced per baryon incorporated into stars. 
Similarly, the ionising efficiency above the atomic cooling mass 
is Ca = /*.a/esc,at7ion,a. The cosmic mass fraction collapsed into 
dark matter haloes above the ionised IGM feedback threshold is 
given by 


7^coll,i (^) 



dMM^iz), 


( 2 ) 


where pc is the critical cosmological density and is the Sheth- 
Tormen mass function. Similarly, the collapsed fractions for atomic 
cooling haloes below Mi and for minihaloes are 


Fco\\,a,(.z') 


7^coll,m (^) 





dMM—iz), 
dn , , 


(3) 

(4) 


We denote the ionised volume of an HII region per unit gas 
mass in a dark matter halo and unit ionising efficiency at redshift « 
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as V{zon, z), where Zau is the redshift corresponding to the forma¬ 
tion of the halo (the total ionising volume for each halo is given by 
V = To determine this value we solve the equation of 

motion of the ionisation front, Ri = (V 


j p3 o 

—f = 3Hiz)Rf + —^-C{z){nH)aBRf, (5) 
at 4:7r [tih} 

where H{z) is the Hubble parameter, as = 2.6 x 
is the case B recombination coefficient of hydrogen at T = 10^ 
K, (n_tr) is the mean cosmic hydrogen density, and C{z) = 
(tt-Hii)/(ttHii)^ is the clumping factor of the ionised IGM. For 
each Mq of star forming gas with ionising efficiency normalised to 
unity, we assume the time-dependent rate of ionising photon emis¬ 
sion, Afy, is 


No {t ^ 10® Vr), 

No X (f/10® ®yr)-4 ® {t > 10® ®yr), 


( 6 ) 


Given these assumptions the LW background is 

, , c(l -I- zf , ,dtH , /, 

JySN\Z) = - — - / dz ), ( 8 ) 

2max 

where c is the speed of light and fn is the Hubble time. The LW 
luminosity per frequency per comoving volume is given by 


eLw{z) = 


V ' 


-??LW,a + 


SFRDn 


rfLW ,m 


(9) 


where SFRD is the star formation rate density, nip is the proton 
mass, rjuw is the number of LW photons per baryon produced in 
stars, i?LW = 1.9 x 10~^^ erg and Ai^lw = 5.8 x 10^^ Hz. 
The subscripts “a” and “m” denote atomic cooling haloes and mini¬ 
haloes as above. We compute SFRDa from the collapsed fractions 
of haloes above the atomic cooling threshold. 


SFRDa(2) = [1 - Q(g)] +Pb/«,a^%^. (10) 


dt 


dt 


where No = 9-25 x 10^^s“^. Over the entire lifetime of the stellar 
population this normalisation yields 1 ionising photon per baryon 
incorporated into stars. The exact form of Aly (t) does not impact 
our results since the majority of the photons are produced much 
faster than the Hubble time for all redshifts relevant to our calcula¬ 
tions. 


2.2 Lyman-Werner feedback 

As described above, the LW background can dissociate molecular 
hydrogen, increasing the minimum mass of minihaloes that host 
Pop III star formation. We incorporate this into our model by self- 
consistently computing the mean LW background, Jlw ( z ), and us¬ 
ing it to set the minimu m mass of minihaloe s hosting Pop III star 
formation. Mm (see, e.g. lHaiman et al. l2000h . We assume the min¬ 
imum mass is equal to 


Mm = 2.5 x 10® 


1 + z 
26 


(l -I- 6.96 (47rJLw(2))® 


(7) 

Hz-1 


where Jlw is in units of 10“ ergs s” 
dFialkov et alj|201^ . This formula gives an increase in minimum 
mass du e to LW radiat i on tha t is consistent with the si mula- 
tions of Machacek et al.l ( 1200 ih . lO’Shea & NormanI ( l2008h . and 
IWise & Abell (l20 07t). The mass for Jlw = 0 is taken as the “op¬ 
timal fit” fro miFialkov et al.] (l2012h which was calibrat ed with the 
simulations of istacv et al. ( 2011 ) and iGreif et al.l (1201 ih . 

We compute Jlw by making a simple “screening” assump¬ 
tion that the IGM is nearly transparent to LW photons until they are 
redshifted into a Lyman series line and absorbed, removing them 
from the LW band. We approximate this by assuming that at a red- 
shift of z all emitted LW photons can be seen from sources out to 
2max = 1.015 X 2 . The factor of 1.015 is used because 1.5 per cent 
is approximately the amount a typical LW photon can be redshifted 
before reaching a Lyman series line. In reality, the LW attenuation 
as a function of frequency will follow a more complicated charac¬ 
teristic “sawtooth” shape. However, we expect our simple approxi¬ 
mation to be reasonably accurate. We find that in our fiducial model 
with constant /*,a described below, our screening reduces the LW 
background by roughly an order of magnitude at 2 « 15 and a fac¬ 
tor of a few at 2 » 30, which is con s istent with more sophisticated 
treatn ients (e.g.lH aiman et ^lTOOd : iRicotti et al.ll200l| ; lAhn et al] 
I 2 OO 9 I: IWolcott-Green et al.ll201 ih. 


Similarly, the minihalo SFRD is given by 

SFRDn,(2) = (1 - Q{Z)) . (11) 

Note that throughout this work, we assume a LW photon escape 
fractio n of unity. This is consistent with the results of ISchauer et al.l 
( l2015h . who find that for a single large Pop III star in a moderate 
size minihalo (6.9 x 10® Mq and 2.1 x 10® Mq ) the “far-field” 
LW escape fraction is unity. For a single star in a halo closer to 
the atomic cooling threshold (1.2 x 10^ Mq), they find an escape 
fraction of « 0.7. Due to the ID nature of these calculations, and 
since we expect a higher rate of star-formation in halos close to 
the atomic cooling threshold, these values should be taken as lower 
limits on the escape fraction. 

To determine the form of Jlw ( 2 ) self-consistently, we itera¬ 
tively compute the entire star formation evolution, reionisation his¬ 
tory, and Jlw( 2 ), using the Jlw( 2 ) computed in previous steps 
until we achieve convergence. 


2.3 Model parameters 

Here we discuss the physical parameters that enter our calculations 
and their fiducial values. Without including LW feedback these pa¬ 
rameters are C(z), Ca, Cm, Mm, Ma, and Mi. With LW feedback 
^ion,a5 ^ion,m5 t?LW,a, and ?7LW,m are also required. 
As explained below, we consider two different models of the star 
formation efficiency in atomic cooling haloes. In the “redshift- 
dependent" model, /*,a varies as a function of cosmic time and 
in the “redshift-independent" model it is constant. We summarize 
the fiducial parameter choices for both of these models in Table 
[T] We treat /*,m as a free parameter and determine how it is con¬ 
strained by Planck in § 3. To demonstrate that these constraints do 
not depend strongly on our fiducial parameter choices we vary each 
of the other parameters subject to the constraint that reionisation is 
completed by 2 = 6 and find that our conclusions remain robust. 

Next we describe the choice for each of our fiducial parame¬ 
ters. We adopt a redshift-dependent clumping factor of the ionised 
IGM parameterised by 

C{z) = 2[^'^ + 1 . ( 12 ) 

This formula is similar to the clumping factor found in the Illus- 
tris simulation for gas below 20 times the mean baryon density 
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4 E. Visbal et al. 


Table 1. Physical parameters used in our reionisation model and their fiducial values. See § 2.3 for more details. 


Parameter 

Description 

Fiducial Value 

C{z) 

clumping parameter 


ea 

ionising efficiency of atomic cooling haloes 

600/*,a [ 400/*,a in redshift-independent model] 

/*,a 

atomic cooling halo star formation efficiency 

see FigureO[0.05 in redshift-independent model] 

Ma 

atomic cooling mass 

5.4 X 10 ’^ 

Mi 

ionised IGM feedback mass 

1.5 X 10* (^)“^-®M© 

Vim, a 

LW photons per stellar baryon in atomic cooling haloes 

4000 

Cm 

ionising efficiency of minihaloes 

40000/*,m 

/*,m 

minihalo star formation efficiency 

treated as free parameter 

Mm 

minimum minihalo mass for star formation via H 2 

2.5 X 10® (1 + 6.96 (47rJLw(z))°'‘‘’’) Mq 

7LW,m 

LW photons per stellar baryon in minihaloes 

80000 



Figure 1. Important mass scales in our reionisation model as a function of 
redshift. For the case with LW feedback, we use the Jlw(^) froni Figure 

a 

dSauer et alJl201.‘I) . lFinlatOT et alJ (l2012h find a clumping factor of 
the ionised IGM similar to this relation as well. 

For the atomic cooling mass, we take a fiducial value of Ma = 
5.4 X 10^ • This corresponds to the typical minimum 

mass of h aloes that are able cool in the absence of H 2 in the simu¬ 
lations of iFernandez et al.] ( l2014h . We set the fiducial ionised IGM 
feedback mass to Mi = 1.5 x 10® ^ This corresponds to 

a halo circular velocity of 20 km s“^, which Ipiikstra et ^ ll2004ll 
find sets the mass scale where feedback becomes important. The 
minimum minihalo mass is assumed to follow Eq.Q which is cali¬ 
brated with the simulations mentioned above. We plot the important 
mass scales discussed here in Figure[T] 

Converting between star formation efficiency and ionisation 
efficiency requires values for 77ic,n and /esc for both atomic cool¬ 
ing haloes and minihaloes. We assume tyion.a = 4000, which cor¬ 
responds to a stellar population with a Salpeter IMF from 0.1 — 
100 Mq and metallicity of Z = 0.0004 (see Table 1 in ISamui et al.l 
(l2007h to see how this quantity changes for different IMFs and 
metallicities). For minihaloes we adopt ryion.m = 80000, which 
is accurate for Pop III stars with masses greater than ~ 200 Mq 
( ISchaerej2002h . Smaller stars would reduce this value (e.g a reduc¬ 
tion by a factor of ~1.40 for 80 Mq stars or ~ 3 for 40 Mq stars). 


For the escape fractions, we take /esc,a = 0.15 in the redshift- 
dependent ft,a. model and /esc,a = 0.1 in the redshift-independent 
ft,a model. We assume /esc,m = 0.5 in minihaloes. These val¬ 
ues ar e generally consistent with the simulations of IWise et al.l 
(l2014h . but we regard this parameter as uncertain. Simulations 
of minihaloes have o ften found escape fra c tions higher than 0. 5 
JWhalen et ^ |2004 iKitavama et ^ |2004 lAlvarez et al.] l2006h . 
suggesting that our fiducial choice is conservative (a low assumed 
value puts weaker limits on Pop III star formation). As discussed 
below, we find that ea(oc /esc,a) cannot be taken to be significantly 
lower than our fiducial model without reionisation occurring at 
z < 6. Higher values would lead to even more stringent constraints 
on ft,in than we present below. When including LW feedback, we 
assume one LW photon per ionising photon (i.e. T^LW.a = 4000 
and ?7LW,m = 80000). This is a reasonably good assu mption for 
a wid e range of IMFs and metallicities (see table 4 in ISchaerej 
I 2 OO 2 I) . Note that only the combinations »7ion,a/esc,a, »7ion,m/esc,m, 
/*,ar?LW,a, and /*,mr?LW,m appear in our model. None of the indi¬ 
vidual parameters appear alone outside of these products. 

We calibrate our fiducial values of Sa and ft,a with observa¬ 
tions of the UV LF at a « 6 and abundance matching. For a given 
absolute UV magnitude, Muv (at a rest-frame wavelength of 1600 
A), we find a corresponding halo mass, M, satisfying 

fM-uy coo 

/ dMuV'/<(Muv) = Eduty / , (13) 

J-00 Jm 


where (j)(M\jy) is th e best fit LF function at 2 « 5.9 from 
ISouwens et alJ bOl.'th (a Schechter function with M, = —20.94, 
(j)t = 0.5 X 10“®, and a = —1.87). We assume that only a frac¬ 
tion of dark matter haloes host bright UV galaxies at a given time. 
This is parameterised with the duty cycle, taken to be Cduty = 0.1 
in the fiducial c ase. This value is consisten t with galaxy clustering 
measurements dBarone-Nugent et alj|2014l) . Once we associate an 
absolute UV magnitude with each halo mass in the relevant range, 
we convert the magnitude to a SFR with the following relation 


SFR 
Mq yr-1 


= 2.24 X 10"^® 


ergs Hz“i 


(14) 


This is the r atio of the dust-correct ed SFRD to the luminosity den¬ 
sity given in lSouwens et al.l ll2015h (see their table 7). This conver¬ 
sion assumes a Salpeter IMF with mass range from 0.1 — 125 Mq 
and solar metallicity. We note that this is not the exact same IMF 
and metallicity used to calibrate r]\oa,a- However, we do not expect 
this to have a large impact on our results and indeed the discrep¬ 
ancy may be justified since the larger observed galaxies are likely 
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Figure 2. SFR per halo mass as a function of halo mass computed with 
abundance matching and the UV LF 5.9 assuming a duty cycle of 10 

per cent. The solid curve shows the masses corresponding to the observed 
data and the dashed line is a power-law extrapolation. The black circle cor¬ 
responds to Ma at 2 = 6. Note that we do not use this extrapolation below 
Ma. 



Figure 3. Star formation efficiency in atomic cooling haloes as a function 
of redshift in our redshift-dependent case computed with Eq. [15] The solid, 
dashed, and dotted curves are for a duty cycle, eduty. of 0.1, 0.5, and 0.02, 
respectively. 

to have somewhat higher metallicity than the smaller more abun¬ 
dant galaxies in atomic cooling haloes driving reionisation. 

In Figure [2| we plot the SFR per halo mass as a function of 
halo mass. The abundance matching calculation indicates that star 
formation is most efficient in a few x 10^° Mq haloes and falls 
off at both lower and higher masses. The observed data extends 
down to a few times 10® M©. Below this value we use a power law 
extrapolation. We emphasise that this approach is conservative, be¬ 
cause it diminishes global star-formation in atomic cooling haloes 
significantly towards higher redshifts, and leaves more room for 
ionising radiation from minihaloes. We use this SFR calibration to 
estimate the star formation efficiency in haloes above the atomic 
cooling threshold, /*,a, by taking the total instantaneous star for¬ 
mation rate divided by the total rate at which gas is being incorpo¬ 


rated into virialised dark matter haloes 

/^‘dMSFR(M)£duty^ 

~ „ df'coll.a • ' 

Pb dt 

Here SFR(M) is the star formation rate in haloes of mass M 
given by our abundance matching calculation. We note that the 
star formation efficiency does not change substantially if we sub¬ 
stitute the limits of integration to correspond to haloes above the 
ionised IGM mass threshold. For the results presented below, we 
consider two different fiducial models of /*,a. In the first, redshift- 
dependent model, we assume that the SFR(M) relationship is 
fixed with respect to redshift and compute how /*,a(2) varies as 
a function of cosmic time from Eq. (see Figure 1^. We note 
that the reduction in /*,a at high redshift is mostly due to the 
evolution of the halo mass function. In particular, the quantity 
( Im' essentially the same redshift 

dependence as /*,a in Figure 0 At high redshifts, the mass in 
haloes relative to the rate of halo collapse is much lower. The ex¬ 
act SFR(M) relation shown in Figure sets the normalization 
of /*,a(2), but only has a small effect on its redshift dependence. 
Figure [3] also shows the dependence of /t,a on Cduty Changing 
the duty cycle from our fiducial value of 0.1 raises or lowers the 
normalization of the calibrated star formation efficiency and has a 
relatively small effect on the redshift dependence. In the redshift- 
independent model we assume /,,a = 0.05, which is the value we 
compute with Eq.[T5]at z = 6. 

2.4 Metal enrichment 

Up to this point, we have made the simplifying assumption that 
minihaloes host only Pop III star formation and atomic cooling 
haloes host only Pop II star formation. However, in reality, super¬ 
novae winds enrich some minihaloes with metals enabling Pop II 
star formation. Conversely, strong LW feedback can lead to the for¬ 
mation of atomic cooling haloes without metals, resulting in Pop III 
star formation. To estimate the impact of these effects, we consider 
a simple idealized model of metal enrichment where the smallest 
star-forming haloes are enriched by Pop III supemovae and subse¬ 
quently form larger enriched haloes. Specifically, we assume that 
Pop III star formation occurs with efficiency /*,iii only between 
Mm and 2Mm, given by Eq.[7] Pop II star formation is assumed 
to form in all larger haloes with efficiency /*,ii. In haloes forming 
Pop II stars we adopt the same parameters (/esc, ?7ion, ?7 lw, etc.) as 
we did for atomic cooling haloes described above (and utilize both 
the redshift-independent and redshift-dependent star formation effi¬ 
ciency models). Similarly, for Pop III star-forming haloes we adopt 
the parameters for minihaloes described above. Operationally, our 
model is computed by changing the upper limit of integration in 
Eq.EI and the lower limit of integration in Eq. to 2Mm ■ When 
LW feedback is included, if the LW background is strong enough 
to suppress star formation in all minihalos (i.e. if the value of Jlw 
causes Mm in Eq. Qto be greater than Ma), we assume that Pop 
III stars form in atomic cooling haloes between Ma and 2Ma (with 
efficiency /*,iii) and Pop II stars form in larger haloes. We model 
this by changing the limits of integration in Eq.|4|to Ma and 2Ma, 
and the lower limit of integration in Eq. |3]to 2Ma. Note that for 
simplicity we have chosen a factor of two in the mass range of ha¬ 
los hosting Pop III stars. This value is not expected to be a precise 
description of the true mass range. However, a factor of two is not 
unreasonable since two halos with masses equal to the minimum 
star forming mass could merge to produce a metal enriched halo. 
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Figure 4. Results of our fiducial reionisation model with the redshift-dependent ft^a,{z) shown in Figure|^ We consider three different models for minihaloes: 
/*,m = 0.001 with (solid curves) and without LW feedback (dashed curves), and a case without any contribution from minihaloes (dotted curves). We plot 
Q{z), SFRD(2) (for minihaloes and atomic cooling haloes), t{z), and ^ for each of these models. In the t{z) (lower left) panel, the thin solid, dashed, 
and dotted lines show the Planck measurement, with the la, and 2a errors, respectively. For the case without LW feedback, early ionisation leads to a high r 
which is inconsistent with Planck. 


Note that this treatment is not self-consistent for very low val¬ 
ues of /»,iii because there will not be enough metals produced in 
haloes smaller than 2M^ to lead to a transition to Pop II stars. 
Thus, for low /,,iii, our model is conservative in the sense that we 
shut-off Pop III star formation in smaller haloes earlier than ex¬ 
pected, leading to weaker constraints on the Pop III star formation 
efficiency. For reference, /»,iii « 10~® in a 10® Mq minihalo cor¬ 
responds to 150 Mq of stars, which could be enough to produce 
a pair instability supernova and enrich gas to the levels required for 
Pop II star formation. 


We also point out that our simple model focuses on “self¬ 
enrichment” of haloes by their progenitor haloes and neglects the 
impact of minihalo enrichment from winds emitted by nearby 
larger haloes. This works in the opposite direction of the inconsis¬ 
tency described above, enriching some small haloes that are metal- 
free in our model. Despite these shortcomings, we expect our sim¬ 
ple treatment to give a rough indication of the impact of metal en¬ 
richment on the reionisation model described above. 


3 RESULTS 

We present the results of our semi-analytic reionisation model in 
the following three subsections. First, we present results from our 
model with the simplifying assumption that minihaloes contain 
only Pop III stars and atomic cooling haloes contain only Pop II 
stars. In Section 3.2, we show how these results could change due 
to metal enrichment in minihaloes with the prescription described 
in Section 2.4. Finally, in Section 3.3 we present limits on the total 
density of Pop III stars formed across cosmic time with and without 
LW feedback and minihalo metal enrichment. 

3.1 No minihalo metal enrichment 

In Figure [T] we show results for our reionisation model with 
redshift-dependent /*,a and the fiducial parameters discussed 
above. We plot Q{z), SFRDa( 2 ), SFRDin(.z), the CMB optical 
depth, t{z), and along with the observational limits from 
Planck. We consider three cases: /*,m = 0.001 without LW feed¬ 
back, /*,in = 0.001 including LW feedback, and /*,in = 0. For the 
case with LW feedback, we plot Jlw(z) in Figure [S] The choice 
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Figure 5. The LW background for our fiducial redshift-dependent /*^a 
model with /*,m = 0.001 in units of 10“^^ ergs s“^ cm“^ Hz“^ Sr“^. 

of star formation efficiency plotted represents moderate formation 
of massive Pop III stars (for /*_m = 0.001, a typical 10® Mq 
minihalo would form 150 Mq in stars). We note that the total 
SFRDa at jg = 6 is sim ilar to that inferred from the observations of 
ISouwens et al.l ^20151) (see their figure 18) since we calibrate with 
their UV LF and their assumed limiting magnitude roughly corre¬ 
sponds to our Mi. However at higher redshifts our SFRD is higher 
somewhat higher because we include the contribution from fainter 
galaxies. 

We compute the optical depth assuming helium is singly 
ionised at the same time as hydrogen and doubly ionised instan¬ 
taneously at z = 3. This leads to 

where ctt is the Thompson scattering cross section, Y = 0.24 and 
X — 0.76 are the helium and hydrogen mass fractions, and r^He = 
1 for z' > 3 and t^hb = 2 for z' ^ 3. 

It is clear from Figure |4] that without LW feedback, a Pop III 
star formation efficiency of = 0.001 in minihaloes is incon¬ 
sistent with the Planck optical depth. Partial reionisation begins 
early, leading to a r more than 3cr higher than the Planck measure¬ 
ment. LW feedback reduces the SFRD in minihaloes by roughly an 
order of magnitude below a = 30. This leads to r = 0.071, which 
is consistent with Planck and leads us to the main conclusion of 
this paper: if massive Pop III stars form in minihaloes with mod¬ 
erate efficiency, LW feedback, metal enrichment, or a comparable 
other suppression of PopIII star-formation is necessary to prevent 
r from being too high compared to Planck data. 

We emphasise that our main conclusions are not sensitive to 
the exact choices of fiducial parameters. To demonstrate this, we 
vary each of the parameters in our model associated with atomic 
cooling haloes, subject to the cons traint that reion isation is essen¬ 
tially complete (Q ~ 1) by 2 « 6 jMesingej|2Q10h . and determine 
how efficiently minihaloes can produce Pop III stars and still be 
consistent with Planck. In Figure|6l we plot r as a function of /,,in 
(assuming 77ion,m = 80000 and/esc,m = 0.5). We plot this both for 
our redshift-independent and redshift-dependent /,,a models with¬ 
out LW feedback and individually vary the other physical param¬ 
eters (besides those associated with minihaloes): ta, Ma, Mi, and 
C{z). We vary these parameters in the direction which lowers r. 


allowing a higher /*,m. Thus, we present conservative upper lim¬ 
its on the efficiency of Pop III stars that can form in minihaloes 
without LW feedback. We also show how r depends on /*,in in the 
fiducial case with LW feedback. Additionally, we plot the redshift 
when reionisation is complete, Zr, for the same models in Figure|7] 
Note that the values plotted for a lower ta are the lowest possible, 
since lowering them any more than this prevents reionisation from 
completing before 2 = 6. 

Overall we find broadly consistent results for the two dif¬ 
ferent parameterisations of our reionisation model. The redshift- 
dependent /*,a model with fiducial parameters and no LW feed¬ 
back requires /*,m ^ 3 x 10“^ to be consistent to within Icr of 
the Planck t . This corresponds to < 50Mq of stars in a typ¬ 
ical 10® M© dark matter halo. Thus, more massive Pop III stars 
forming in most minihaloes without LW feedback is inconsistent 
with Planck. We find that changing the fiducial parameters does not 
have a large effect on these results. Even for the variations in C{z), 
which we regard as extremely conservative, the limits on /*,in only 
change by about a factor of 2. When LW feedback is included, r 
is consistent with Planck for /*,m < 0.002. This would permit 
massive Pop III stars in typical minihaloes, before LW feedback 
prevents star formation. 

We find similar, but even more severe constraints on /*,m in 
the redshift-independent /*,a model without LW feedback. How¬ 
ever, when LW feedback is included, because the background is 
higher at early times due to the greater /*,a, the star formation 
efficiency in minihaloes can be /*,in < 0.003 without violating 
the Planck constraints. Note that we do not show the case with 
C{z) = 6 (^4^) + 1 because this leads to a reionisation history 

where reionisation is not complete until after 2 = 6. 

3.2 Impact of minihalo metal enrichment 

In Figure]^ we plot r and 2 r as a function of Pop III star formation 
efficiency, /*,iii, for our reionisation model including the treatment 
of metal enrichment described in Section 2.4 with and without LW 
feedback. We find that for both the redshift-dependent and redshift- 
independent /*,ii models, when LW feedback and metal enrich¬ 
ment are included, the Planck l-cr limit on r corresponds to a Pop 
III star formation efficiency of /»,iii ~ 10~®. For the redshift- 
independent /,,ii case with LW feedback and metal enrichment, 
the limit on Pop III star formation efficiency goes down by a factor 
of ~ 4 compared to the case with LW feedback and no metal en¬ 
richment. This is because our treatment of metal enrichment causes 
Pop III stars to form in the smallest atomic cooling haloes, increas¬ 
ing the total amount of ionizing photons produced for sufficiently 
high Pop III star formation efficiency. 

For the redshift-dependent /*,ii case with LW feedback, we 
find that the limits on Pop III star formation efficiency are similar 
with or without metal enrichment. With ~ 10“®, we find that 
in both cases, for 2 < 20, the LW feedback sets the minimum mini¬ 
halo mass to be roughly half the atomic cooling mass, which leads 
to Pop III star formation in the same mass range with or without 
metal enrichment and explains the close similarity in r. In general, 
we find that the effects of metal enrichment plus LW feedback do 
not greatly change the constraints on Pop III star formation effi¬ 
ciency obtained with LW feedback alone. 

We note that in reality metal enrichment will not lead to a 
global change in the mass range of halos that host Pop III stars. In¬ 
stead there will be a complex interplay between radiative feedback 
and metal enrichment leading to the mass range varying strongly 
as a function of position. As such, we caution the reader that the 
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Figure 6. The optical depth due to electron scattering vs. minihalo star-formation efficiency for two different parameterisations of our reionisation model, 
and for various different parameter values. For comparison, we plot the r measured by Planck (solid black line) and the Icr and 2cr error bars (dashed and 
dotted black lines). For all cases we have assumed /esc ,m = 0.5 and = 80000. The redshift-dependent /*_a model (see Figure|3) is plotted in the 

left panel and the redshift-independent, /*,a = 0.05, model is plotted in the right panel. For the cases without LW feedback all parameters have been varied 
in the direction which decreases r, providing a conservative upper limit on /*,m. The variations in ta represent the largest possible decrease without causing 
reionisation to complete at 2 : < 6 (in the right panel, we do not include the factor of 3 change in C (2 = 6) because this delays reionisation to 2 < 6). 




Figure 7. Same as Figure|6] but showing the redshift when reionisation is completed, 2 r, instead of r. The black circles denote the star formation efficiency 
conesponding to the l-cr Planck optical depth limits. 


results in this subsection are meant only to give a rough indication 
of the possible impact of metal enrichment on the constraints on 
Pop III stars. 


3.3 Limits on total density of Pop III stars 

Next, we examine how the Planck optical depth limits the total den¬ 
sity of Pop III stars formed over cosmic time. We compute this by 
integrating SFRDui{z) (SFRDiii{z) in the case with metal en¬ 
richment) with respect to cosmic time. 

In Figure we plot the cumulative density of Pop III stars, 
p*.iii, with /*, ,m (/*, III in cases with metal enrichment) corre¬ 
sponding to the 1-a limits on the Planck optical depth. We show re¬ 
sults for the redshift-independent and redshift-dependent /*,a mod¬ 
els with and without LW feedback and our treatment of metal en¬ 


richment. Interestingly, we find that including LW feedback and/or 
metal enrichment does not have a large impact on the limits of the 
total number of Pop III stars produced. For the redshift-dependent 
and redshift-independent /*_a cases, we find total limiting Pop III 
densities of ~ 10® Mq Mpc“® and ~ 3 x 10"* Mq Mpc“®, re¬ 
spectively. Including or not including LW feedback and metal en¬ 
richment does not change this result by more than a factor of a few. 
LW feedback or prompt metal-enrichment reduce the number of 
mini-halos that can form Pop III stars, allowing the efficiency of 
Pop III star formation to rise in those halos, but ultimately the strict 
new limits on r only permit a small number of Pop III stars to form 
overall. 

We note that in our models with LW feedback, hut no 
metal enrichment, Pop III star formation can he completely 
suppressed when the LW background is sufficiency high. This 
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Figure 8 . The CMB optical depth and redshift of reionisation for our semi-analytic model with metal enrichment. The horizontal solid, dashed, and dotted 
lines show the Planck t measurement, l-cr, and 2-cr limits, respectively. The black circles in the right panel indicate the star formation efficiency comesponding 
to the l-CT Planck limits. We do not include the redshift-independent /* jj case without LW feedback, because unphysical Pop II star formation occurs at very 
high redshift in minihaloes leading to an artificially high optical depth. 


can be seen in the right panel of Figure^ where Pop III star for¬ 
mation is halted at z « 15 (however it resumes at 2 ; « 10 when 
the LW background decreases). For the other models, Pop III 
star formation is only completely stopped by the completion of 
reionization. 


4 LIMITS ON POP III BH SEED GROWTH 

If massive Pop III stars are produced in minihaloes, some will end 
their lives as black holes. As these black holes grow, they are ex¬ 
pected to produce X-rays which provide an additional source of 
ionising radiation not included in the model described above. In 
this section, we perform a simple calculation to estimate the im¬ 
pact of black hole accretion on the IGM. 

To estimate the total mass of black holes produced, we use the 
fiducial redshift-dependent /*,a model described above, including 
LW feedback and /*,m = 0.001. We assume that 10 per cent of 
the stellar mass formed in minihaloes ends up as black holes (com¬ 
puted with the SFRD in Eq.[TTJ- This is approximately the fraction 
obtained for a Salpeter IMF with mass limits of 1 Mq and 100 Mq 
assuming that stars with with initial mass between 40 Mq and 100 
Mq collapse directly to black holes. Note that a more top heavy 
IMF could produce a somewhat larger fraction. 

We assume that black holes grow at the Eddington limit a frac¬ 
tion EBH.duty of the time with radiative efficiency Cr = 0.1. For the 
entire population of black holes this results in an accretion rate 

= 2.2 x10-'’£bh, duty i^PBHyr-', (17) 

where pbh is the total comoving black hole density. If the black 
hole density reaches pbh = 10® MqMpc”®, we turn off accre- 
tion by hand (mim icking inefficient growth or a self-regulation; 
iTanaka et alj|2012n to prevent the density from greatly exceeding 
that of SBMHs in the local universe. 

We compute the effect of the black hole growth on the IGM 
by taking the number of ionisations per volume per time as 

driion _ fioji CrC^ dpBH 

dt ^ E-, dT’ ' ’ 


where E-y = 13.6 eV is energy required to ionise a hydrogen atom. 
An X-ray produced through black hole accretion will ionise a hy¬ 
drogen or helium atom producing a high energy electron. A frac¬ 
tion of this electron’s energy, /ion, will go into producing additional 
ionisations. The value of /ion depends on the energy of the electron 
and the ionised fracti on of the IGM. We e s timate /ion by inter¬ 
polating the results of iFurlanetto & Stoeved j201(t) and assuming 
a typical electron energy of 1 keV. For an ionised fraction close 
to zero this gives us /ion ~ 0.4 and is reduced to nearly zero as 
the ionised fraction approaches unity. Using this rate of ionising 
photon production and the fiducial C{z) described above, we com¬ 
pute the ionised fraction from black hole accretion alone, Qy,{z). 
We assume that because of the large mean free path of X-rays, the 
IGM is unifor mly ionis ed, as opposed to having an ionised bubble 
topology (e.g. lOhir200lll . Since we are only performing a rough es¬ 
timation of the impact of black hole accretion, we do not attempt to 
self-consistently model the combination of X-rays and UV photons 
from stars simultaneously. Note that we also do not model possi¬ 
ble self-reg ulation of black hole growth via X-ray feedback, as was 
explored in iTanaka et alj ( l2012h . 

In Figure [TOl we plot pm{z) and Q^{z) for esH.duty = 
1, 0.1 and 0.01. It is clear that for a duty cycle of 1 the ionisa¬ 
tion is inconsistent with Planck. The optical depth from Qx alone is 
r = 0.122, and taking the total ionisation as the sum Q + Qx (with 
a max of 1) yields r = 0.15. For eBH,duty = 0.1 and 0.01, taking 
the total ionisation as Q + Qx gives r = 0.091 and r = 0.075, 
respectively. Thus, only for esH.duty 0.01 is this model compat¬ 
ible with the la Planck limits. 


5 DISCUSSION AND CONCLUSIONS 

We have performed analytic calculations of reionisation including 
separate contributions from Pop II stars in atomic cooling haloes 
and Pop III stars in minihaloes with and without LW feedback. For 
haloes above the atomic cooling threshold we considered two dif¬ 
ferent models, one with a constant star formation efficiency and the 
other with a constant SFR-halo mass relation leading to a redshift- 
dependent star formation efficiency. We calibrate both of these 
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Figure 9. The total cumulative density of Pop III stars as a function of redshift for Pop III star formation efficiency corresponding to the 1 -cr Planck upper 
limits, including LW feedback and/or metal enrichment. The left (right) panel shows the results for the redshift-dependent (redshift-independent) /,,a model. 


models with the observed z « 6 UV LF. We also incorporated a 
simple treatment of metal enrichment in minihaloes due to Pop III 
supernovae. 

Without LW feedback or metal enrichment, a minihalo star 
formation efficiency greater than ~a few x 10 ““^ creates an early 
partial reionisation incompatible with the Planck optical depth (as¬ 
suming /esc,m = 0.5). In a 10®M 0 minihalo, this star formation 
efficiency corresponds to 50Mq of stellar mass. Thus, typical 
minihaloes could not form more massive Pop III stars. To demon¬ 
strate the robustness of this conclusion, we vary each model pa¬ 
rameter, aside from those associated with minihaloes, and find that 
our conclusions do not change significantly. We note that the ion¬ 
ising efficiency of atomic cooling haloes cannot be lowered signif¬ 
icantly from the fiducial value or reionisation will not be complete 
by z = 6 , inconsistent with measurements of the Lya forest. With 
our empirical calibration of /*,a and 77ion,a = 4000, the escape 
fraction cannot be less than /esc,a ~ 0.1 or reionisation will occur 
too late. 

When we include our self-consistent treatment of LW feed¬ 
back and/or metal enrichment we find that the star formation effi¬ 
ciency in (the massive, M > Mm) minihaloes can be much higher 
without violating the Planck constraints: ft,in ~ a few x 10~®. 
This leads us to the main conclusion of our paper. Without LW 
feedback as strong as the model used in our paper or metal en¬ 
richment, massive Pop III stars in minihaloes will lead to an op¬ 
tical depth of the CMB that is inconsistent with Planck observa¬ 
tions. Another important conclusion of our paper is that the limit 
on the total density of large Pop III stars formed over cosmic 
time due to the 1-a Planck constraints is roughly in the range 
of 10"^ — 10® Mq Mpc“®, irrespective of the feedback prescrip¬ 
tion (i.e. LW, LW-l-metals, etc). For reference, we note that this is 
1.6 X 10“® — 1.6 X 10~® of the total baryon density. We also point 
out that this fraction does not correspond to the density of Pop III 
remnants remaining today because most of the stars have lifetimes 
much shorter than the age of the Universe. 

For our cases without LW feedback or metal enrichment al¬ 
though we quote constraints in terms of /,,m, the relevant quantity 
is Cm = /*,m?7ion,m/esc,m- Thus, if the escape fraction were lower 
than our assumed /esc,m = 0.5 by some factor, the corresponding 
limit on /*,m would go up by the same factor. 


Note that in our model we have ignored the bary on-dark mat¬ 
ter streaming velocity jTseliakhovich & Hiratal2oT^ . At high red- 
shifts this effect can reduce the efficiency of star formation in mini¬ 
haloes. However, since most of the contribution to r in our model 
comes from 2 < 20 , this effect would onl y reduce the star for ma- 
tion efficiency by a factor of a few at most jFialkov et alj|2012h . 

We have also considered how X-rays emitted due to the ac¬ 
cretion of black hole remnants from Pop III stars would impact 
the IGM. We performed a simple estimate of X-ray ionisation 
from black holes produced by our fiducial model with /*,in = 
0.001 and found that unless the duty cycle of black hole ac¬ 
cretion is eBH,duty 0 . 01 , early ionisation produces an opti¬ 

cal depth greater than the Planck Icr limits. While we empha¬ 
sise that our rough estimate is somewhat model dependent (e.g. 
a very hard X-ray spectrum could lead to free-streaming of X- 
rays and weaker constraints), the result is intriguing. To form the 
first super massive black holes (SMBHs), which are more mas¬ 
sive than 10® Mq at 2 « 6, would require that Pop III remnants 
grow at the Eddington limit with a duty cycle of nearly unity. 
This suggests that either some type of feedback (e.g. lTanaka et al.l 
l 20 l^ acts on most, but not all Pop III remnants if they are 
the seeds of the first SMBHs or that SMBHs are seeded by 
a different mechanism such as direct collapse black holes (e.g. 
[Brornm & Loeb||2003|: |yblonteri & Begelma nl20ld : |Diikstra et al.l 
20l4IVisbal et alj2014l:|lnavoshi et alj2015h. 
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Figure 10. The comoving black hole density from remnants of Pop III stars, the fraction of the IGM ionised by X-rays emitted due to black hole accretion, and 
the corresponding t and ^. Note that the ionisation includes no contribution from stars, while r and ^ assume a total ionisation of Q -|- Qx (with a max 
value of 1). The details of the calculation are explained in § 4. The solid, dashed, and dotted curves denote eBH,duty = !■ 0.1, and 0.001, respectively. The 
horizontal solid, dashed, and dotted lines in the t panel denote the Planck measurement, and the Icr and 2 <t errors. Black hole growth in the eeH.duty = 1 
case is stopped by hand at 2 ^ 11 in order to avoid overproducing Pbh- This causes the IGM to begin recombining at 2 < 11 - while somewhat ad-hoc, this 
is conservative since it limits the BH’s contribution to r (see text). 
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